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Abstract 



Inclusive distributions of charged particles in hadronic W decays are experi- 
mentally investigated using the statistics collected by the DELPHI experiment 
at LEP during 1997, 1998 and 1999, at centre-of-mass energies from 183 to 
around 200 GeV. The possible effects of interconnection between the hadronic 
decays of two Ws are not observed. 

Measurements of the average multiplicity for charged and identified particles in 
qq and WW events at centre-of-mass energies from 130 to 200 GeV and in W 
decays are presented. 

The results on the average multiplicity of identified particles and on the posi- 
tion £* of the maximum of the £ p = — log(%) distribution are compared with 
predictions of JETSET and MLLA calculations. 



(Eur. Phys. J. C18(2000)203; Erratum ibid. C25(2002)493) 
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1 Introduction 

The study of the production properties of charged and identified hadrons (7r + , K + , 
K°, p and A) 1 in qq events at LEP 2 allows QCD based models to be tested through the 
comparison with LEP 1 data, in particular the validity of high-energy extrapolations of 
Monte Carlo models tuned at the Z°. 

In the case of WW events in which both W bosons decay hadronically, these kind of 
studies, besides providing checks for QCD-inspired models, are expected to give insights 
into possible correlations and/or final state interactions between the decay products of 
the two W bosons. 



Hadron production in e + e and QCD 

The way quarks and gluons transform into hadrons is complex and can not be com- 
pletely described by QCD. In the Monte Carlo simulations (as in [1]), the hadronisation 
of a qq pair is split into 3 phases. In a first phase, gluon emission and parton branching 
of the original qq pair take place. It is believed that this phase can be described by 
perturbative QCD (most of the calculations have been performed in leading logarithmic 
approximation LLA). In a second phase, at a certain virtuality cut-off scale Q , where 
a s (Q ) is still small, quarks and gluons produced in the first phase are clustered in colour 
singlets and transform into mesons and baryons. Only phenomenological models, which 
need to be tuned to the data, are available to describe this stage of fragmentation; the 
models most frequently used in e + e~ annihilations are based on string and cluster frag- 
mentation. In the third phase, the unstable states decay into hadrons which can be 
observed and identified in the detector. These models account correctly for most of the 
features of the qq events such as, for instance, the average multiplicity and inclusive 
momentum spectra. 

A different and purely analytical approach (see e.g. [2] and references therein) giv- 
ing quantitative predictions of hadronic spectra are QCD calculations using the so-called 
Modified Leading Logarithmic Approximation (MLLA) under the assumption of Local 
Parton Hadron Duality (LPHD) [3,4]. In this picture multi-hadron production is de- 
scribed by a parton cascade, and the virtuality cut-off Qo is lowered to values of the 
order of 100 MeV, comparable to the hadron masses; it is assumed that the results ob- 
tained for partons are proportional to the corresponding quantities for hadrons. 

The MLLA/LPHD predictions involve three parameters: an effective scale parame- 
ter A e fj, a virtuality cut-off Qq in the evolution of the parton cascade and an overall 
normalisation factor -Klphd- The momentum spectra of hadrons can be calculated as 
functions of the variable £ p = — ln(-^=), with p being the particle's momentum and y/s 
the centre-of-mass energy: 

^-=K LP}1D .f^ p ,X,X) (1) 

with 

A- = log^ ; A = log* (2) 

Due to uncertainties from higher order corrections A e g- cannot be identified with A^^. 
In equation (1), n is the average multiplicity per bin of £ p , and the function / has the 
form of a "hump-backed plateau". It can be approximated by a distorted Gaussian [5,6] 

1 Unless otherwise stated antiparticlcs arc included as well. 
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£>G(£; N, £, a, s k , k) = -j= exp Qfc + i Sfc 5 - 1(2 + A;)5 2 + l Sfc 5 3 + ^k5 4 ^ , (3) 

where 5 = (£ — £)/o", £ is the mean of the distribution, a is the square root of its variance, 
5^ its skewness and k its kurtosis. For an ordinary Gaussian these last two parameters 
vanish. The mean, £, coincides with the peak of the distribution, £*, only up to next-to- 
leading order in a s . 

To check the validity of the MLLA/LPHD approach, one can study the evolution of 
the position of the maximum, £*, as a function of ^/s. In the context of MLLA/LPHD 
the dependence of £* on the centre-of-mass energy can be expressed as [2,5]: 



C = y {^ + VcjY - C/Y^j + F h 



(A), (4) 



where 



v^/2 / lliV c /3 + 2n / /(3^) \ 2 / N c \ 

F " l0g XT ' °~{ W c ) -\UN e /3-2n f /3j> (5) 

with N c being the number of colours and iif the number of active quark flavours in 
the fragmentation process. The function Fh(X) depends on the hadron type through 
A = log(<5o/A e fr) [2], and it can be approximated as F h (X) = — 1.46A + 0.207A 2 with an 
error of ±0.06. 



Interference and final state interactions in W decays 

The possible presence of interference due to colour reconnection and Bose-Einstein 
correlations (see for example [7-12] and [13,14] for reviews) in hadronic decays of WW 
pairs may provide information on hadron formation at a time scale smaller than 1 fm/c. 
At the same time it can induce a systematic uncertainty on the W mass measurement in 
the 4-jet mode [13] comparable with the expected accuracy of the measurement at LEP 2. 

Interconnection can happen due to the fact that the lifetime of the W (tw — fr/^w — 
0.1 fm/c) is an order of magnitude smaller than the typical hadronisation times. The 
interconnection between the products of the hadronic decays of different Ws in WW pair 
events can occur at several stages: (1) from colour rearrangement between the quarks 
coming from the primary branching, (2) due to gluon exchanges during the parton cas- 
cade, (3) in the mixing of identical pions or kaons due to Bose-Einstein correlations. The 
first two are QCD effects. They can mix the two colour singlets and produce hadrons 
which cannot be uniquely assigned to either W. The perturbative effects are colour sup- 
pressed and the possible shift is expected to be only about 5 MeV in the W mass [7]. 

Non-perturbative effects need model calculations. Several models have been proposed 
(for reviews see for example [13] and [15]) and have already been included in the widely 
used event generators PYTHIA [1], ARIADNE [16] and HERWIG [17]. In these models 
the final state quarks after the parton shower can be rearranged to form colour singlets 
with probabilities which in some cases are free parameters. The shift on the W mass in 
these models is typically smaller than 50 MeV [18], but other observables are affected 
by colour rearrangement. Generally these models suggest a small effect on the total 
charged particle multiplicity, of the order of —1% to —2% [15,19,14]. Dedicated detailed 
simulations of the response of the DELPHI detector to such models showed that this effect 
is substantially unaffected by the event selection criteria and by the detector performance. 
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For identified heavy particles, such as K + and p , the effects due to colour reconnection 
are expected to be stronger [20] , but the experimental verification is complicated by losses 
in statistics. The same applies to the particle spectrum at low momentum [15]. Bose- 
Einstein Correlations could also slightly change the multiplicity for (4g) events in some 
models [12,21]. 

The WW events allow a comparison of the characteristics of the W hadronic decays 
when both Ws decay in hadronic modes (referred to here as the (4g) mode) with the 
case in which only one W decays hadronically ((2q) mode). These characteristics should 
be the same in the absence of interference between the hadronic decay products from 
different W bosons. 

Previous experimental results based on the statistics collected by LEP experiments 
during 1997 (see [22-24] for reviews) did not indicate at that level of statistics the presence 
of interconnection or correlation effects. 

This paper presents measurements of: 

• the charged particle multiplicities for the qq events in the data sample collected by 
the DELPHI experiment at LEP during 1997, 1998 and 1999, at the centre-of-mass 
energies from 183, 189 and from 192 to 200 GeV respectively; 

• the charged particle multiplicity and inclusive distributions for WW events at 183 
and 189 GeV (multiplicity values at 189 GeV are expected to be slightly higher than 
that at 183 GeV due to increased phase space. However this effect is below the 
precision obtainable with the present data samples); 

• The average multiplicities for identified charged and neutral particles K + , K°, 
p and A), and the position £* of the maximum of the £ p distribution for identified 
particles in qq events from 130 GeV to 189 GeV and in WW events at 189 GeV. 

2 Data Sample and Event Preselection at 183 and 
189 GeV 

Data corresponding to total luminosities of 157.7 pb -1 (54.1 pb -1 ) at centre-of-mass 
energies around 189 (183) GeV , and data taken in 1999 corresponding to total luminosi- 
ties of 25.8 pb -1 , 77.4 pb -1 and 83.8 pb -1 at centre-of-mass energies around 192 GeV, 
196 GeV and 200 GeV respectively were analysed. A description of the DELPHI detector 
can be found in reference [25]; its performance is discussed in reference [26]. 

A preselection of hadronic events was made, requiring at least 6 charged particles and 
a total transverse energy of all the particles above 20% of the centre-of-mass energy y/s. 
In the calculation of the energies E, all charged particles were assumed to have the pion 
mass. Charged particles were required to have momentum p above 100 MeV/c and below 
1.5 times the beam energy, a relative error on the momentum measurement Ap/p < 1, 
angle 9 with respect to the beam direction between 20° and 160°, and a distance of closest 
approach to the interaction point less than 4 cm in the plane perpendicular to the beam 
axis (2 cm in the analyses of identified charged particles) and less than 4/sin# cm along 
the beam axis (2 cm in the analyses of identified charged particles). 

After the event selection charged particles were also required to have a track length 
of at least 30 cm, and in the charged identified particles analysis a momentum p > 200 
MeV/c. 

The influence of the detector on the analysis was studied with the full DELPHI simula- 
tion program, DELSIM [26]; events were generated with PYTHIA 5.7, using the JETSET 
fragmentation with Parton Shower (PS) [1] with parameters tuned to fit LEP 1 data from 
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DELPHI [27]. The initial state for the WW 1999 sample was generated using EXCAL- 
IBUR version 1.08 [28]. The particles were followed through the detailed geometry of 
DELPHI with simulated digitizations in each detector. These data were processed with 
the same reconstruction and analysis programs as the real data. 

To check the ability of the simulation to model the efficiency for the reconstruction 
of charged particles, the samples collected at the Z° pole during 1998 and 1997 were 
used. From these samples, by integrating the distribution of £ E = — m (^|)> where E 
is the energy of the particle, corrected bin by bin using the simulation, the average 
charged particle multiplicities at the Z° were measured. The values were found to be 
20.93 ± 0.03(stat) and 20.60 ± 0.03(stat) respectively, in satisfactory agreement with the 
world average of 21.00 ± 0.13 [29]. The ratios of the world average value to the measured 
multiplicities at the Z°, respectively 1.0033 ±0.0064 and 1.0194 ±0.0065 from the Z° data 
in 1998 and 1997, were used to correct the measured multiplicities at high energies in the 
respective years. 

The cross-section for e + e~ — > qq(7) above the Z° peak is dominated by radiative qq7 
events; the initial state radiated photons (ISR photons) are generally aligned along the 
beam direction and not detected. In order to compute the hadronic centre-of-mass energy, 
y/s*, the procedure described in reference [30] was used. The procedure clusters the 
particles into jets using the DURHAM algorithm [31], excluding candidate ISR photons 
and using a y cut = 0.002. The reconstructed jets and additional ISR photons are then 
fitted with a three constraint fit (energy and transverse momentum, leaving free the z 
component of the missing momentum). The hadronic centre-of-mass energy, y/s', is the 
invariant mass of the jets using the fitted jet energies and directions. 

3 Analysis of charged particles in qq events 

3.1 Centre-of-mass energies of 183 and 189 GeV 

Events with ^fs' / y/s above 0.9 were used to compute the multiplicities. A total of 3444 
(1297) hadronic events were selected from the data at 189 (183) GeV, by requiring that 
the multiplicity for charged particles was larger than 9, that the total transverse energy 
of the charged particles exceeded 0.2a/s, and that the narrow jet broadening [32] was 
smaller than 0.065. From the simulation it was calculated that the expected background 
coming from WW and Z°Z° decays was 432±52 (127±21) events. The contamination 
from double radiative returns to the Z°, within 10 GeV of the nominal Z° mass, was 
estimated by simulation to be below 5%. Other contaminations (from Z°ee, Wez/, 77 
interactions and Bhabhas) are below 2% in total. 

The average multiplicity of charged particles with p > 0.1 GeV/c measured in the 
selected events at 189 GeV (183 GeV), after subtraction of the WW and Z°Z° backgrounds 
estimated by simulation, was 24.58 ± 0.16(stat) (23.96 ± 0.23(stat)), to be compared to 
24.52 ±0.05(stat) (24.30 ± 0.07 (stat)) in the qq PS simulation including detector effects. 
The dispersion (square root of the variance) of the multiplicity distribution in the data 
was 7.57 ± O.ll(stat) (7.00 ± 0.16(stat)), to be compared to the dispersion from the qq 
PS simulation of 7.24 ± 0.03 (7.20 ± 0.05{stat)). 

Detector effects and selection biases were corrected for using a qq simulation from 
PYTHIA with the JETSET fragmentation tuned by DELPHI without inital state radia- 
tion. The corrected average charge multiplicity was found to be < n >= 27.37±0.18(stat) 
(< n >= 26.56 ± 0.26(stat)), and the dispersion was found to be D = 8.77 ± 0.13(stat) 
(D = 8.08 ± 0.19(stat)). 
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The average multiplicity was computed by integrating the distribution, since the 
detection efficiency depends mostly on the momentum of the particle, after correcting for 
detector effects bin by bin using the simulation. The £g distribution was integrated up 
to a value of 6.3, and the extrapolation to the region above this cut was based on the 
simulation at the generator level. 

After multiplying by the Z° corrections factors from section 2, the following values 
were obtained: 



These values include the products of the decays of particles with lifetime r < 10 9 s. 
The systematic errors were obtained by adding in quadrature: 

1. the propagated uncertainty of the average values in the Z° correction factors, ±0.18 
(±0.17) for the multiplicity. 

2. the effect of the cuts for the reduction of the background. The value of the cut 
on the narrow jet broadening was varied from 0.045 to 0.085 in steps of 0.010, in 
order to estimate the systematic error associated with the procedure of removing 
the contribution from WW events. The new values for the average charged particle 
multiplicity and the dispersion were stable within these variations, and half of the 
difference between the extreme values, 0.07 and 0.07 (0.06 and 0.12) respectively, 
were added in quadrature to the systematic error. The effect of the uncertainty on 
the WW cross-section was found to be negligible. 

3. the uncertainty on the modelling of the detector response in the forward region. The 
analysis was repeated by varying the polar angle acceptance of charged particles 
from 10-170 degrees to 40-140 degrees, both in the high energy samples and in the 
computation of the Z° correction factors. The spread of the different values obtained 
for the multiplicities and for the dispersions were found to be respectively 0.18 and 
0.08 (0.23 and 0.03). The effect of the variation of other track selection criteria was 
found to be negligible, and the same applies to the higher centre-of-mass energies. 

4. the systematic errors due to the statistics of the simulated samples, 0.04 (0.06) for 
the multiplicity and 0.04 (0.06) for the dispersion. 

5. the uncertainty on the calculation of the efficiency correction factors in the multi- 
plicity. The values of the multiplicities, before applying the Z° correction factors, 
were also estimated: 

• from the observed multiplicity distribution as 27.37 (26.56); 

• from the integral of the rapidity distribution (with respect to the thrust axis), 

Ut = \ In g^ 11 is the absolute value of the momentum component on the 

thrust axis) as 27.43 (26.59). 

Half of the differences between the maximum and the minimum values of the mul- 
tiplicity calculated from the multiplicity distribution and from the integration of 
the ut and £e distributions, 0.03 in both cases, were added in quadrature to the 
systematic error. 

6. Half of the extrapolated multiplicity in the high-^ region, 0.14 (0.12). 



(n)is9GcV = 27.47 ± 0.18(stat) ± 0.30(syst) 
£>i89 Gcv = 8.77 ± 0.13(stat) ± O.ll(syst) 
MissGeV = 27.05 ± 0.27(stat) ± 0.32{syst) 
£>i83 Gcv = 8.08 ± 0.19(staf) ± 0.U(syst) . 



(6) 
(7) 
(8) 
(9) 
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Table 1: Number of events selected, number of events expected for the signal (qq7) 
and for the background (WW and ZZ), estimated from simulation, and the measured 
multiplicities and dispersions for the three different energies. 

As a cross-check, a simulated sample based on HERWIG plus DELSIM was also used 
to unfold the data; the results were consistent with those based on PYTHIA plus DELSIM 
within the statistical error associated to the size of the Monte Carlo sample. 

3.2 Centre-of-mass energies of 192 to 200 GeV 

To check the ability of the simulation to model the efficiency for the reconstruction 
of charged particles, the sample collected at the Z° calibration runs of 1999 was used, 
following the procedure described in section 2. The average charged particle multiplicity 
at the Z° was measured to be 20.82 ± 0.03(stat), in satisfactory agreement with the 
world average. The ratio of the world average value to the measured multiplicity at the 
Z°, 1.0084 ± 0.0064, was used to correct the measured multiplicities at centre-of-mass 
energies of 192 to 200 GeV. 

For each of the energies a separate analysis was performed following the procedure 
described in the previous subsection. 

The number of events selected, the number of expected signal and background events, 
estimated with Monte Carlo simulation, and the measured multiplicities and dispersions 
are listed in Table 1 for the centre-of-mass energies of 192 to 200 GeV. The systematic 
errors were estimated as in 3.1; a breakdown is shown in the table (the numbering of the 
sources of systematic error corresponds to the one in the previous subsection). 

These results were then combined at an average centre-of-mass energy of 200 GeV, 
according to the following procedure. First each result was rescaled to a centre-of-mass 
energy of 200 GeV, calculating the scaling factors from the simulation. Then a weighted 
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average was computed using the statistical error as a weight. The systematic error is 
taken as the weighted average of the systematic errors, increased (in quadrature) by the 
difference between the values obtained when rescaling and when not rescaling to 200 GeV. 
This gives: 

(n) 20 0GcV = 27.58 ± O.lQ(stat) ± 0A5(syst) (fO) 
£>2O0GeV = 8.64 ± O.lS(stat) ± 0.20{syst) . (11) 

As a cross-check, a simulated sample based on HERWIG plus DELSIM was also used 
to unfold the data; the results were consistent with those based on PYTH1A plus DELSIM 
within the statistical error associated to the size of the Monte Carlo sample. 

4 Classification of the WW Events and Charged 
Multiplicity Measurement 

About 4/9 of the WW events are WW — > qiq^qsqi events. At threshold, their topology 
is that of two pairs of back-to-back jets, with no missing energy; the constrained invariant 
mass of two jet-jet systems is close to the W mass. Even at 183 and 189 GeV these 
characteristics allow a clean selection. 

Another 4/9 of the WW events are WW — > qiq~2 tv events. At threshold, their topology 
is 2-jets back-to-back, with a lepton and missing energy opposite to it; the constrained 
invariant mass of the jet-jet system and of the lepton-missing energy system equals the 
W mass. 



4.1 Fully Hadronic Channel (WW — > qiq2q3q4) 

Events with both Ws decaying into qq are characterised by high multiplicity, large 
visible energy, and tendency of the particles to be grouped in 4 jets. The background is 
dominated by qq(7) events. 

The events were pre-selected by requiring at least 12 charged particles (with p > 
100 MeV/c), with a total transverse energy (charged plus neutral) above 20% of the 
centre-of-mass energy. To remove the radiative hadronic events, the effective hadronic 
centre-of-mass energy computed as described in section 2, was required to be above 
110 GeV. 

The particles in the event were then clustered to 4 jets using the LUCLUS algorithm [1], 
and the events were kept if all jets had multiplicity (charged plus neutral) larger than 
3. It was also required that the separation between the jets (enjoin value) be larger than 
6 GeV/ c. The combination of these two cuts removed most of the remaining semi-leptonic 
WW decays and the 2-jet and 3-jet events of the qqy background. 

A five constraint fit was applied, imposing energy and momentum conservation and 
the equality of two di-jet masses. Of the three fits obtained by permutation of the jets, 
the one with the smallest \ 2 was selected. Events were accepted only if 

d sc1 = f""*™" > Q Q04rad QeV _ 1 

^ 'max 'max -^min) 

where E min and E max are respectively the smallest and the largest fitted jet energy, and 
m i n is the smallest angle between the fitted jet directions. The details of the selection 
variable D se \ can be found in [33]. The purity and the efficiency of the selected data 
sample from the 189 (183) GeV data were estimated using simulation to be about 76% 
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and 80% (75% and 80%) respectively. The data sample consists of 1256 (427) events, 
where 1255 (422) were expected from the simulation. The expected background was 
subtracted bin by bin from the observed distributions, which were then corrected bin by 
bin using scaling factors computed from the simulation generated using PYTHIA with 
the JETSET fragmentation tuned by DELPHI (EXCALIBUR plus JETSET for the 1999 
data) without initial state radiation. 

Finally, the average multiplicity of charged particles (n^) was estimated by inte- 
grating the £e distribution up to a value of 6.3 (and estimated above this value with 
simulation) and multiplying by the Z° correction factors from section 2. The following 
values were obtained: 



The systematic errors account for: 

1. The propagated uncertainty of the average values in the Z° correction factors, ±0.24 
(±0.24). 

2. The spread of the measured values from the reference values by redoing the analysis 
varying the selection criteria, 0.05 (0.11). 

3. Modelling of the detector in the forward region. The analysis was repeated by 
varying the polar angle acceptance of charged particles from 10-170 degrees to 40- 
140 degrees, both in the WW samples and in the computation of the Z° correction 
factors. The spreads of the different measured values were found to be 0.01 (0.08). 

4. Limited statistics in the simulated sample 0.03 (0.05). 

5. Variation of the qqy cross-sections within 5%: 0.01 (0.01). 

6. Calculation of the correction factors. The value of (n^), before applying the Z° 
correction factors, was also estimated: 

• from the observed multiplicity distribution as 38.96 (37.39); 

• from the integral of the rapidity distribution (with respect to the thrust axis) 
as 39.16 (36.97). 

• from the integral of the pt distribution (with respect to the thrust axis) as 38.95 



Half of the difference between the maximum and the minimum value, 0.11 (0.23), 
was added in quadrature to the systematic error. 
7. Uncertainty on the modelling of the 4-jets qq background, 0.05 (0.14). The uncer- 
tainty on the modelling of this background is the sum in quadrature of two contri- 
butions: 

• Uncertainty on the modelling of the 4-jet rate. The agreement between data 
and simulation was studied in a sample of 4-jet events at the Z°, selected with 
the DURHAM algorithm for y cut ranging from 0.003 to 0.005. The rate of 4-jet 
events in the simulated sample was found to reproduce the data within 10%. 
The correction due to background subtraction was correspondingly varied by 
10%, which gives an uncertainty of 0.00 (0.01). 

• Uncertainty on the multiplicity in 4-jet events. The average multiplicity in 4-jet 
events selected at the Z° data in 1998 (1997), with the DURHAM algorithm 
for a value of y cut = 0.005, is larger by 0.9% ± 0.3%(stat) (2.0% ± 0.5%(stat)) 
than the corresponding value in the simulation. A shift by 0.9% (2.0%) in 
the multiplicity for 4-jet events induces a shift of 0.05 (0.14) on the value in 
Equation (12) (Equation (13)). 



(n (4<?) )i89Gev = 39.12 ± 0.33(stat) ±0.36(syst) 
(™ (4<?) >i83GeV = 38.11 ± 0.57(stat) ± O.U(syst) . 



(12) 
(13) 



(37.42). 
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8. Half of the extrapolated multiplicity in the high-^g region, 0.23 (0.20). 

The presence of interference between the jets coming from the different Ws could 
create subtle effects, such as to make the application of the fit imposing equal masses 
inadequate. For this reason a different four constraint fit was performed, leaving the 
di-jet masses free and imposing energy-momentum conservation. Of the three possible 
combinations of the four jets into WW pairs, the one with minimum mass difference was 
selected. No x 2 cu t was imposed in this case. The average multiplicity obtained was 
again fully consistent (within the statistical error) with the one measured in the standard 
analysis. 

A simulated sample based on HERWIG plus DELSIM was also used to unfold the 
data; the results were consistent with those based on PYTHIA plus DELSIM within the 
statistical error associated to the size of the Monte Carlo sample. 

The distribution of the observed charged particle multiplicity in (Aq) is shown in 
Figure lc(la). 

The value of the corrected multiplicity in the low momentum range 0.1 to 1. GeV/c, 
where the interconnection effects are expected to be most important, was found to be 
14.47±0.20 (13.67±0.34), where the errors are statistical only. 

After correcting for detector effects, the dispersion was found to be: 



In the systematic error: 

1. 0.10 (0.15) accounts for the spreads of the measured values from the reference value 
when varying the event selection criteria; 

2. 0.03 (0.01) is due to the modelling of the detector in the forward region. The 
dispersions were also measured using only charged particles with polar angle between 
40 and 140 degrees and the differences with respect to the reference value were 
considered in the systematic error; 

3. 0.03 (0.05) from the limited simulation statistics. 

4.2 Mixed Hadronic and Leptonic Final States (WW— ► qiq^f) 

Events in which one W decays into lepton plus neutrino and the other one into quark 
and antiquark are characterised by two hadronic jets, one energetic isolated charged 
lepton, and missing momentum resulting from the neutrino. The main backgrounds to 
these events are radiative qq production and four-fermion final states containing two 
quarks and two oppositely charged leptons of the same flavour. 

Events were selected by requiring seven or more charged particles, with a total energy 
(charged plus neutral) above 0.2^ and a missing momentum larger than 0.1^/s. Events 
in the qqy final state with ISR photons at small polar angles, which would be lost inside 
the beam pipe, were suppressed by requiring the polar angle of the missing momentum 
vector to satisfy | cos^^l < 0.94. 

Including the missing momentum as an additional massless neutral particle (the can- 
didate neutrino), the particles in the event were clustered to 4 jets using the DURHAM 
algorithm. The jet for which the fractional jet energy carried by the highest momentum 
charged particle was greatest was considered as the "lepton jet" . The most energetic 
charged particle in the lepton jet was taken as the lepton candidate, and the event was 




8.72 ± 0.23(stat) ± 0.11(syst) 
8.53 ± 0.39{stat) ± 0.16{syst) . 



(14) 
(15) 
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rejected if its momentum was smaller than 10 GeV/c. The neutrino was taken to cor- 
respond to the missing momentum. The event was discarded if the invariant mass of 
the event (excluding the lepton candidates) was smaller than 20 GeV/c 2 or larger than 
110 GeV/c 2 . 

At this point three alternative topologies were considered: 

• Muon sample: if the lepton candidate was tagged as a muon and its isolation angle, 
with respect to other charged particles above 1 GeV/c, was above 10°, the event 
was accepted either if the lepton momentum was greater than 20 GeV/c, or if it 
was greater than 10 GeV/c and the value of the parameter required by the 
DURHAM algorithm to force the event from a 3-jet to a 4-jet configuration was 
greater than 0.003. 

• Electron sample: if the lepton candidate had associated electromagnetic energy de- 
posited in the calorimeters larger than 20 GeV, and an isolation angle (defined as 
above) greater than 10°, the event was accepted if the required value of yf!^ 4 was 
greater than 0.003. 

• Inclusive sample: the events were also accepted if the momentum of the lepton 
candidate was larger than 20 GeV/c and greater than half of its associated energy, 
the missing momentum was larger than 0.1y/s, the required value of was greater 
than 0.003, and no other charged particle above 1 GeV/c existed in the lepton jet. 

The purity and the efficiency of the selected data sample from the 189 (183) GeV data 
were estimated using simulation to be about 89% and 54% (88% and 55%) respectively. 
The data sample consists of 633 (256) events, where 689 (235) were expected from the 
simulation. The expected background was subtracted bin by bin from the observed dis- 
tributions, which were then corrected bin by bin using scaling factors computed from the 
simulation generated using PYTHIA with the JETSET fragmentation tuned by DELPHI 
(EXCALIBUR plus JETSET for the 1999 data) without initial state radiation. 

By integrating the £e distribution up to a value of 6.3 (and estimated above this 
value from simulation at generator level), and multiplying by the Z° correction factors 
from section 2, the following values were obtained for the charged multiplicity for one W 
decaying hadronically in a WW event with mixed hadronic and leptonic final states: 

(n {2q) )is9GcV = 19.49 ± OM(stat) ±0.27(syst) (16) 
(n (2,?) )i83GeV = 19.78 ± 0A9(stat) ± 0A3(syst) . (17) 

In the systematic error: 

1. 0.12 (0.12) accounts for the propagated uncertainty of the average values in the Z° 
correction factors; 

2. 0.16 (0.16) accounts for the spreads from the reference values when changing the 
event selection criteria; 

3. 0.05 (0.02) accounts for modelling of the detector in the forward region. The analysis 
was repeated by varying the polar angle acceptance of charged particles from 10-170 
degrees to 40-140 degrees, both in the WW samples and in the computation of the 
Z° correction factors. The spreads of the different measured values were found to be 
0.05 (0.02); 

4. 0.03 (0.05) accounts for limited statistics in the simulated samples; 

5. the variation of the qq(7) cross-sections within 5%, gives a negligible contribution 
to the systematic error; 

6. 0.12 (0.36) accounts for the uncertainty on the correction factors. The value of 
(n^), before applying the Z° correction factors, was also estimated: 
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• from the observed multiplicity distribution as 19.53 (19.60); 

• from the integral of the rapidity distribution (with respect to the thrust axis) 
as 19.30 (18.88). 

• from the integral of the pt distribution (with respect to the thrust axis) as 19.39 
(19.36). 

Half of the difference between the maximum and the minimum value, 0.12 (0.36), 
was included in the systematic error. 
7. Half of the extrapolated multiplicity in the high-^ s region, 0.12 (0.10). 

A simulated sample based on HERWIG plus DELSIM was also used to unfold the 
data; the results were consistent with those based on PYTHIA plus DELSIM within the 
statistical error associated to the size of the Monte Carlo sample. 

The distribution of the observed charged particle multiplicity in (2q) is shown in 
Figure ld(lb). 

The value of the corrected multiplicity in the low momentum range 0.1 to 1. GeV/c 
was found to be 7.29±0.19 (7.15±0.28) (where the errors are statistical only). 
After correcting for detector effects, the dispersions were found to be: 

D^J GeY = 6A9±0.21(stat)±0A3(syst) (18) 
£>gJ GeV = 6.51 ± 0.33(stat) ± 0.25(syst) . (19) 

In the systematic error: 

1. 0.14 (0.08) accounts for the variation of the cuts; 

2. 0.41 (0.23) accounts for the modelling of the detector in the forward region. The 
dispersions were also measured using only charged particles with polar angle between 
40 and 140 degrees and the differences with respect to the reference value were 
considered in the systematic error; 

3. 0.03 (0.05) is due to the limited simulation statistics. 



5 Analysis of Interconnection Effects from Charged 
Particle Multiplicity and Inclusive Distributions 

Most models predict that, in case of colour reconnection, the ratio between the mul- 
tiplicity in (4g) events and twice the multiplicity in (2q) events would be smaller than 1; 
the difference is expected to be at the percent level. It was measured: 

V2(n( 2 <?)) ) 

( (g) \ 
\2(n( 2 i)) ) 

In the calculation of the systematic error on the ratio, the correlations between the 
sources of systematic error were taken into account. If the systematic errors are taken as 
uncorrelated, except for the errors on the Z° correction factors and the modelling of the 
detector in the forward region, for which full correlation is assumed, a compatible value 
of ±0.014 (±0.022) is obtained for the systematic error. 



= 1.004 ± 0M8(stat) ± O.OU(syst) (20) 
= 0.963 ± 0.028(stat) ± 0M5(syst) . (21) 
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Using as weights the inverse of the sum in quadrature of the statistical and systematic 
errors, one obtains a weighted average of 

/ n (4<?)\ \ 

N (2 ' = 0.990 ± 0M5(stat) ± OMl(syst). (22) 

In the presence of interconnection, the deficit of multiplicity is expected to be con- 
centrated in the region of low momentum. The corrected momentum distributions in 
the (4g) and in the (2g) cases are shown in Figures 2 and 3 (Figures 4 and 5 show the 
distributions in terms of the £g variable). The systematic error in the momentum region 
between 0.1 and 1 GeV/c was explicitely recomputed. We measure: 

/ ( 4 g)\ 0.1<p<lGeV/c 

= 0.992 ± 0.029(stat) ±0MQ(syst) (23) 

189 GeV 
0.1<p<l GeV/c 

= 0.956 ± 0.044(stat) ±0.022{syst) . (24) 

183 GeV 



2(n( 2 <?)) 



2(n( 2 9)> 



Using as weights the inverse of the sum in quadrature of the statistical and systematic 
errors, one obtains a weighted average of 

( n (4q)\ 0.1<p<lGeV/c 

= 0.981 ±0.024(stat) ±0M3(syst). (25) 



Other inclusive distributions (rapidity and p T for example) taking into account the 
orientation of particles could display a larger sensitivity with respect to interconnection 
effects. Special care should be taken, since the definition of the thrust axis could introduce 
a bias between {2q) and (4q) events. To estimate the effect of this bias, the following 
procedure was used. First the distributions of rapidity and pr were computed as in the 
case of the momentum for the (2q) and the (4g). Then, a set of (4g)-like events was 
constructed by mixing pairs of (2q) events (constructed (4q) sample). To construct this 
sample the hadronically decaying Ws in (2q) events were boosted back to their rest frame. 
For each real (4g) event, each W was then replaced by a W from a (2g), boosted forward 
in such a way that the directions of the jet axis were the same. 

If there would be no bias from the definition of the thrust axis (or if the correction 
for detector effects could correctly account for the bias), there would be no difference 
between the yx and pt distributions in the constructed (4g) sample and twice the (2g) 
sample. This difference has thus been taken as an estimator of the systematic error from 
the bias, and added in quadrature to the (4g) distribution, to twice the (2q) distribution, 
and to their difference. The distributions of transverse momentum with respect to the 
thrust axis are shown in Figures 6 and 7. It can be seen that the difference between the 
(4g) distribution and twice the (2q) at 183 GeV is concentrated in the low-py region; 
however, part of the difference is due to the fact that the (2q) distribution at low pt lies 
above the simulation. The rapidity distribution, instead, does not display any particular 
feature. 

The dispersion in (4g) events is consistent at both energies with \/2 times the dispersion 
in (2g) events (in the calculation of the systematic error on the ratio, the correlations 
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between the sources of systematic error were taken into account): 

= 0.95 ± 0.04(stat) ± 0.07(syst) 

/ £)(4?) \ 

-= = 0.93 ± 0.06(stat) ± 0.04(syst) . 



Using as weights the inverse of the sum in quadrature of the statistical and systematic 
errors, one obtains a weighted average of 

/ £)(4«) \ 

\ V2DW ) = °' 94 ± °-° 3(stat) ± °- 03 (^ st )- ( 26 ) 

In conclusion, no depletion of the multiplicity was observed in fully hadronic WW 
events with respect to twice the semileptonic events, at this statistical precision; a possible 
depletion at the percent level can however not be excluded. 



6 Identified Particles from e + e — > qq 

This section describes the results obtained for n + , K + , K°, p and A with data recorded 
by DELPHI at LEP 2. After the description of the event selection for identified particles 
at energies up to 189 GeV, the additional criteria for hadron identification are described. 



6.1 Event selection at 130 and 136 GeV 

After the hadronic preselection described in Section 2, events with \fs' > 0.85a/s were 
used for further analysis. Data recorded at these two energies were combined and are 
referred to as the 133 GeV sample. The total of ~ 12 pb -1 recorded by DELPHI yields 
1387 events while 1406 are expected from simulation. 



6.2 Event selection at 161 and 172 GeV 

Selected events at 161 GeV were required to have a minimum of 8 and a maximum of 
40 charged particles, \fs' > 0.85a/s and a visible energy of at least 50% of y/s. A cut was 
imposed on the polar angle 9 of the thrust axis to select events well within the acceptance 
of the detector. It was found that a selection based on the narrow jet broadening, B Nl is 
effective in removing the WW events and minimises the bias introduced on the remaining 
event sample. At threshold about 30 WW events are expected. Selecting events with 
.B/v < 0.12 reduces this background by 50%. 

Using this selection 342 events are expected from simulation, while 357 were selected 
from the data, with an estimated remaining WW background of 15 events. 

At 172 GeV, in addition to the criteria described above, events were required to have 
at most 38 charged particles and Bn < 0.1. This leads to 267 selected events, with 264 
expected from simulation, out of which 36 are WW background. 



6.3 Event selection at 183 and 189 GeV 



The event selection for charged identified particles at 189 (183) GeV follows very closely 
the procedure already described in sections 2 and 3. Events with ^fs' / y/s above 0.9 and 
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more than 9(8) charged particles with p > 200 MeV/c were used. WW background was 
suppressed by demanding B N < 0.1(0.08). A total of 3617(1122) events were selected 
with an expected background from WW and Z°Z° of 789 (146) events. 

6.4 Selection of charged particles for identification 

A further selection was applied to the charged particle sample to obtain well identi- 
fiable particles. Two different momentum regions were considered, above and below 0.7 
GeV/c, which correspond respectively to the separation of samples identified solely by 
the ionization loss in the Time Projection Chamber (TPC) and by the Cerenkov detec- 
tors (RICH) and TPC together. Below 0.7 GeV/c tighter cuts were applied, namely to 
eliminate secondary protons. There had to be at least 30 wire hits in the TPC associated 
with the track and the measured track length had be to larger than 100 cm. In addition 
it was required at least two associated VD layer hits in r<p and an impact parameter in 
the rcj) plane of less then 0.1 cm. If there were less than two associated VD layers in 
z, the corresponding impact parameter had to be less than 1 cm, else less then 0.1 cm. 
Particles above 0.7 GeV/c were required to have a measured track length bigger than 30 
cm and good RICH quality, i.e. presence of primary ionization in the veto regions. Only 
particles which were well contained in the barrel region of DELPHI (| cos(#) |< 0.7) were 
accepted. 

6.5 Analysis 

For an efficient identification of charged particles over the full momentum region, 
information from the ionization loss in the TPC ("dE/dx") and information from the 
DELPHI RICH detectors were combined, using dedicated software packages [26]. One 
package fine-tunes the Monte Carlo simulation concerning detector related effects (such 
as slight fluctuations in pressures and refractive indices, background arising from photon 
feedback, crosstalk between the MWPC readout strips, <5-rays, track ionization photo- 
electrons, etc.), and another package derives identification likelihoods from the specific 
energy loss, the number of reconstructed photons and the mean reconstructed Cerenkov 
angles respectively. The likelihoods are then multiplied and rescaled to one. From these, 
a set of "tags" which indicate the likelihood for a particular mass hypothesis (ir, K, and 
p) are derived. Throughout this analysis leptons were not separated from pions. Their 
contribution to the pion sample was subtracted using simulation. 

A matrix inversion formalism was used to calculate the true particle rates in the 
detector from the tagged rates. The 3x3 efficiency matrix is defined by 



where type i, j can be either of n,K ± ,p(p). It establishes the connection between the true 
particles in the RICH/TPC and the tagged ones: 




Number of type i hadrons tagged as type j hadrons 



(27) 



Number of type i hadrons 




(28) 



The inverse of the efficiency matrix works on the three sets of tagged particles in two 
ways. First a particle can have multiple tags, meaning that the information from the 
tagging is ambiguous. This is not unlikely because in this analysis the low statistics of 
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the data samples force rather loose selection criteria to be applied. Secondly a particle 
can escape identification. Both effects can be corrected by this method. The average 
identification efficiency is approximately 85% for pions and 60% for kaons and protons, 
whereas the purities are approximately 85% and 60% respectively. They show a strong 
momentum dependence. 

K° and A candidates were reconstructed by their decay in flight into 7r + 7r~ and p7r~ 
respectively. Secondary decays candidates, V°, in the selected sample of hadronic events 
were found by considering all pairs of oppositely charged particles. The vertex defined 
by each pair was determined so that the x 2 °f the hypothesis of a common vertex was 
minimised. The particles were then refitted to the common vertex. The selection criteria 
were the "standard" ones described in [26]. The average detection efficiency from this 
procedure is about 36% for the decay K° — > it + it~ and about 28% for the decay A — > 
in multi-hadronic events. The background under the invariant mass peaks was subtracted 
separately for each bin of V° momentum. The background was estimated from the data 
by linearly interpolating two sidebands in invariant mass: 

• between 0.40 and 0.45 GeV/c 2 and between 0.55 and 0.60 GeV/c 2 for the K° ; 

• between 1.08 and 1.10 GeV/c 2 and between 1.14 and 1.18 GeV/c 2 for the A. 

6.6 Calibration of the efficiency matrix using Z° data 

The Z° data recorded during each year for calibration were used to tune the above 
described matrix before applying it to high energy data. This is made possible by the 
fact that studies at the Z° pole [34] established that the exclusive particle spectra are 
reproduced to a very high level of accuracy by the DELPHI-tuned version of the genera- 
tors. Therefore deviations of the rates of tagged particles between data and simulation in 
this sample can be interpreted as detector effects. Comparison of the tag rates allows a 
validation of the efficiency matrix, which would be impossible to measure from the data 
due to the limited LEP 2 statistics. The matrix is corrected so that it reproduces the 
simulated rates, assuming that the correction factors are linear in the number of tagged 
hadrons. The discrepancies were found to be smaller than 4%. This is taken into account 
when calculating the systematic uncertainties. 

As the high energy events are recorded over a long time period, stability of the identi- 
fication devices becomes a major concern. Variations in the refractive index or the drift 
velocity in the RICH detectors may significantly change the performance of the identifica- 
tion. To estimate the effect of these variations on the measurement, the radiative returns 
to the Z° were used. Such events were selected among the events passing the hadronic 
preselection, by requiring in addition that they contained at least 8 charged particles, 
they had \fs' < 130 GeV, and a total energy transverse to the beam axis of more than 
30 GeV. Figure 8 shows the good agreement for differential cross-sections for this event 
sample which may be taken as an indication of the stability of the detector during the 
year at the few percent level. 

6.7 £ p distributions and average multiplicities 

After background subtraction, the tagged particle fractions were unfolded using the 
calibrated matrix. The full covariance matrix was calculated for the tag rates using 
multinomial statistics. It was then propagated to the true rates of identified particles 
using the unfolding matrix. 



1G 



Vs [GeV] 


Particle 


(n) 






PYTHIA 5.7 


HERWIG 5.8 


ARIADNE 4.8 


Data 


133 


7r± 
K± 
K° 

PP 
A 


19.90 

2.37 

2.40 

1.11 

0.34 


19.97 

2.32 
2.64 
0.89 
0.49 


19.64 
2.30 

1.29 


19.84 ±0.29 ±0.44 
2.60 ±0.26 ±0.13 
2.51 ± 0.21 ± 0.11 
1.56 ±0.25 ±0.09 
0.50 ±0.07 ±0.03 


161 


7r± 
K± 
K° 

PP 
A 


21.24 
2.63 
2.56 
1.26 


21.76 
2.44 
2.78 
0.96 


21.12 
2.44 

1.40 


20.75 ±0.58 ±0.50 
2.87 ±0.55 ±0.25 
2.65 ± 0.34 ± 0.14 
1.21 ±0.48 ±0.02 


172 


K± 
K° 

PP 
A 


21.77 
2.68 


22.31 
2.48 


21.68 
2.50 


21.79 ±0.68 ±0.47 
2.09 ±0.74 ±0.29 




1.30 


0.99 


1.45 


1.78 ±0.73 ±0.25 


183 


K± 
K° 

PP 
A 


22.28 
2.74 
2.66 
1.33 
0.39 


22.82 
2.53 
2.91 
1.00 
0.55 


21.18 
2.55 

1.48 


21.79 ±0.36 ±0.46 
2.83 ±0.37 ±0.13 
1.81 ± 0.14 ± 0.10 
1.32 ±0.34 ±0.17 
0.33 ±0.04 ±0.03 


189 


K± 
K° 

PP 
A 


22.56 

2.77 

2.67 

1.35 

0.39 


23.10 
2.55 
2.93 
1.02 
0.56 


22.47 
2.57 

1.50 


22.19 ±0.24 ±0.46 
3.15 ±0.21 ±0.24 
2.10 ±0.12 ±0.10 
1.19 ±0.17 ±0.41 
0.40 ±0.03 ±0.03 



Table 2: Average multiplicities of particles 7r + , K + , K°, p and A at 133, 161, 172, 183 
and 189 GeV. The first uncertainty is statistical and the second is systematic. 

The £ p distribution was corrected bin by bin for detector acceptance and selection 
efficiency, using the full detector simulation. As an example, the corrected £ p distributions 
for charged pions, charged kaons and protons at 189 GeV are shown in Figures 9,10,11 
for qq, (4q) and (2g) events respectively. In the figures the predictions from the DELPHI 
simulation as well as a fit to expression (3), if the fit converged, are also shown. Figure 12 
shows the same distributions for K° and A in qq events at 183 and 189 GeV. Within the 
statistics of the data samples analysed, the shapes of the £ p distributions are reasonably 
well described by the generators, with the exception of the K° for which the agreement 
is poorer. 

The multiplicity of the identified final states per hadronic event was obtained by 
integration of the corresponding £ p distributions. The results are shown in table 2. The 
numbers given for charged identified hadrons include decay products from particles with 
a lifetime r < 10~ 9 s. These numbers are compared with the predictions from PYTHIA 
5.7, HERWIG 5.8 and ARIADNE 4.8. 

The following sources of systematic uncertainties are taken into account. 

1. Uncertainties due to particle identification. 



17 



They are mainly due to the uncertainties in the modelling of the detector response. In 
addition for the charged particles there are time dependent effects such as variations 
of the drift velocity in the RICH detectors. 

The unfolding matrix was adjusted using the Z°-calibration data recorded in the 
beginning of data taking in 1996, 1997 and 1998, as well as from the peak period in 
1995. This is described in section 6.6. For 1997 and 1998 data also radiative return 
Z°-events were used. These are better suited as they were recorded under the same 
conditions as the signal events. 

Spectra were obtained using the original and the adjusted matrices. The difference 
between the results obtained was taken as the corresponding systematic uncertainty. 
The relative size of this uncertainty, averaged over the centre-of-mass energies is 
0.006, 0.076 and 0.136 for charged pions, kaons and protons respectively. 
A relative systematic uncertainty of 3% was estimated for K° [35] and 5% for A [36]. 

2. Size of the subtracted WW background (where applicable). 

Variation of the selection criteria results in a 10% uncertainty on the fraction of 
the W contamination in the qq sample. This corresponds to a 5% uncertainty in 
the W cross-section and the size of the W background has been varied accordingly. 
The maximal variation observed in the distributions has been taken as the corre- 
sponding systematic uncertainty. The relative size of this uncertainty, averaged over 
the centre-of-mass energies is 0.003, 0.010 and 0.003 for charged pions, kaons and 
protons respectively, and is 0.014 (0.028, negligible) and 0.032 (0.057) for K° and A, 
respectively, at centre of mass energies of 189 (183, 161) GeV. 

3. Particles with momenta below 0.2 GeV/c or above 50 GeV/c were not identified. 
Their contribution was extrapolated from the simulation. Half of the extrapolated 
multiplicity was added in quadrature to the systematic uncertainty. The relative 
size of this uncertainty, averaged over the centre-of-mass energies is 0.021, 0.003 and 
0.001 for charged pions, kaons and protons respectively, while for K° and A at 133 
GeV, 183 GeV and 189 GeV, is 0.034, and for K° at 161 GeV is 0.042. 

4. For the pions this analysis relies on a subtraction of the lepton contamination using 
simulation. 

An extra uncertainty of 10% of the simulation prediction for the total number of 
leptons is added. The relative size of this uncertainty, averaged over the centre-of- 
mass energies is 0.003. 

In Figure 13 the results are compared to the predictions from PYTHIA 5.7 and HER- 
WIG 5.8. The results shown for energies below 133 GeV (open squares) were extracted 
from reference [29]. The models account for the measurements, with the exception of the 
HERWIG predictions for K° and A and of the PYTHIA predictions for K° at high energy 
which are above the measured values. 

6.8 and its evolution 

An interesting aspect of the ^-distribution is the evolution of its peak position £* 
with increasing centre-of-mass energy. It is determined by fitting a parametrisation of 
the distribution to the peak region. 

One such parametrisation is the distorted Gaussian in equation (3). Another 
parametrisation is a standard Gaussian distribution. While being a more crude approxi- 
mation, it facilitates the analysis in the case of limited statistics. 

Since equation (3) is expected to describe well only the peak region, the fit range has 
to be carefully chosen around the peak. Table 3 shows the results and the fit range, 
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with the x 2 P er degree of freedom for the gaussian fit, except for the K° and A at 133 
GeV where only the distorted gaussian fit was performed. The errors in the data are the 
sum in quadrature of the statistical and systematic uncertainties (generator values were 
extracted from samples of one million events generated with the DELPHI tuned versions 
of the programs) . 

The systematic uncertainty has the following contributions which were added in 
quadrature to the statistical uncertainty of the fit. 

1. Uncertainty of the background evaluation (above the W threshold). 

This source was evaluated as the maximal difference obtained by a variation of the 
WW background cross-section. 

The relative size of this uncertainty, averaged over the centre-of-mass energies is 
0.005, 0.010, 0.002, 0.024 and 0.023 for charged pions, kaons, protons, K° and A 
respectively. 

2. Uncertainty due to the particle identification. 

The analysis was repeated using the calibrated matrices or changes by 3% (K°) or 
5% (A) in the bin contents of the £ p distribution, and the fit redone. The maximal 
differences in the position of the peak thus obtained were added in quadrature. 
The relative size of this uncertainty, averaged over the centre-of-mass energies is 
0.011, 0.017, 0.002, 0.019 and 0.017 for charged pions, kaons, protons, K° and A 
respectively. 

3. Stability of the fit and dependence on the fit range. 

To estimate this effect, which arises from the combination of the limited statistics, the 
resulting need to choose a coarse binning, and the choice of the fit range, systematic 
shifts have been imposed on the data by variation within the statistical uncertainty. 
One standard deviation has been added to the values left of the peak and one 
standard deviation has been subtracted from the values to its right and vice versa. 
The maximum variation is taken as the contribution to the systematic uncertainty. 
The relative size of this uncertainty, averaged over the centre-of-mass energies is 
0.026, 0.112, 0.161, 0.163 and 0.147 for charged pions, kaons, protons, K° and A 
respectively. 

Figure 14 shows the £* values from the Gaussian fits as a function of the centre-of- 
mass energy. The data up to centre-of-mass energies of 91 GeV were taken from previous 
measurements [37]. The fits to expression (4), with Y and C = 0.351 defined as in 
equation (5) and Fh and A e fj taken as free parameters, were done separately for each 
particle type and are superimposed on the data points (lines). Figure 14 shows that 
(within the statistics of the data samples analysed) the fitted functions follow well the 
data points. This suggests that MLLA+LPHD gives a good description of the observed 
particle spectra. From table 3 it is shown that there is fair agreement between the data 
and the predictions from the generators (PYTHIA 5.7, HERWIG 5.8, and ARIADNE 4.8). 

7 Identified hadrons in WW events 

The selections of (4g) and (2q) WW events in the analysis of identified hadrons in 
WW events at 189 GeV are similar to the procedures described in [38]. A feed forward 
neural network, trained with back-propagation on PYTHIA simulated events, was used 
to improve the separation of (4g) WW events from 2-fermion (mainly Z°/7 — ■> qq(7)) and 
4-fermion background (mainly Z°Z° events). The network input variables, the training 







fit range 


PYTHIA 


HERWIG 


ARIADNE 


Data 






Gauss. dist. G. 


Gauss, dist. G. 


Gauss, dist. G. 


Gauss, dist. G. 


Gauss. X^/ndf dist. G. 


133 


K± 
K° 

PP 
A 


[3.3 : 4.9] [2.5 : 5.5] 
[2.1 : 3.3] 

[0.6 : 5.4] 

[2.1 : 3.3] 

[0.6 : 4.8] 


3.90 4.11 
2.85 2.99 

2.87 
3.04 3.09 

2.79 


3.90 4.12 
2.86 2.99 

3.15 
3.04 3.10 

2.99 


4.00 4.15 
2.85 2.99 

3.04 3.04 


4.04 ±0.13 0.017 4.32 ±0.37 
2.90 ±0.30 1.39 

1.63 2.86 ±0.43 
2.74 ±0.33 1.86 

3.22 2.81 ±0.66 


161 


tt* 
K± 
K° 

PP 
A 


[3.5 : 5.1] [2.6 : 4.5] 
[2.1 : 3.7] 

[2.1 : 3.7] 


4.09 4.17 
2.99 2.98 

3.01 3.14 


4.03 4.08 
3.00 3.00 

3.20 


4.12 4.18 
3.01 2.98 

3.14 3.06 


4.11 ±0.15 1.10 4.14 ±0.48 
3.16 ±0.50 0.57 

3.08 ±0.84 1.1 


172 


7r± 
K± 
K° 

PP 
A 


[3.5 : 5.1] [2.6 : 5.5] 
[2.1 : 3.7] 

[2.1 : 3.7] 


4.13 4.20 
3.04 3.02 

3.16 3.14 


4.06 4.12 
3.03 3.01 

3.23 3.41 


4.16 4.21 
3.05 3.01 

3.18 3.12 


4.11 ±0.14 1.07 4.45 ±0.40 
2.98 ±0.84 1.76 

3.55 ±1.10 0.98 


183 


7r± 
K± 
K° 

PP 
A 


[2.6 : 4.5] [3.5 : 5.1] 
[2.1 : 3.7] 

[1.8 : 4.8] [1.2 : 5.4] 
[2.3 : 3.7] [2.1 : 3.7] 
[1.8 : 4.8] [0.6 : 6.6] 


4.16 4.43 

3.08 3.18 

3.09 3.19 
3.21 3.32 
3.02 3.01 


4.10 4.30 
3.06 3.15 

3.31 


4.19 4.44 
3.10 3.23 

3.21 3.39 


4.23 ±0.13 2.8 4.63 ±0.20 
3.08 ±0.29 1.19 
3.33 ±0.63 2.75 3.35 ± 0.63 
3.29 ±0.18 0.98 3.29 ± 0.52 
3.25 ±0.51 0.87 3.22 ± 0.51 


189 


K± 
K° 

PP 
A 


[2.6 : 4.5] 
[2.1 : 4.7] 

[1.8 : 4.8] [0.6 : 6.0] 
[2.1 : 4.7] 

[1.8 : 4.8] [0.6 : 4.8] 


4.22 
3.17 

3.13 3.02 
3.11 

3.02 2.98 


4.08 
3.08 

3.20 




4.18 ±0.10 2.32 

3.10 ±0.25 0.26 

3.06 ±0.48 1.66 2.81 ± 0.43 

3.18 ±0.12 0.01 

3.18 ±0.48 0.38 3.21 ± 0.49 
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fit range 


WW simulation 


data 


X 2 /n.d.f. 




[2.7 : 4.5] 


3.66 


3.65 ±0.02 ±0.24 


1.20/2 




[1.10 : 3.86] 


2.70 


2.61 ±0.09 ±0.33 


0.02 


PP 


[1.10 : 3.86] 


2.73 


2.86 ±0.11 ±0.44 


0.08 



Table 4: Values for £* for 7r ± , and protons from a Gaussian fit to the semi-leptonic 
W data at 189 GeV. The fit was made to the spectra in the W rest frame. The first error 
is statistical and the second reflects the total systematic uncertainties. Also indicated is 
the x 2 1 /n.d.f '. 

procedure and its performance are described in detail in [38]. The purities and the 
selection efficiencies were determined from the simulation to be respectively 73% and 
83% in the (4q) sample and 89% and 75% in the (2q) sample. 

7.1 £ p distributions for identified particles at Mw 

The mass of the W boson is the relevant energy scale which enters in the £ p variable; 
it is thus appropriate to compute the £ p distributions and their peak positions from W 
decays after a boost to the rest frame of the W. The boosting procedure relies on the 
simulation for the part of the spectrum which was not accessible before the boosting, and 
thus introduces a certain model dependence. This does not affect significantly the peak 
region, which is used to extract £*. Figure 11 shows the £ p distributions obtained from 
semi-leptonic W events. 

Table 4 shows the values of £* for the semi-leptonic W data at 189 GeV and for the 
WW simulation. The £* values were obtained by a Gaussian fit in the indicated region. 
The systematic uncertainties were estimated in the same way as described in the previous 
section. The relative uncertainty due to particle identification is 0.03, 0.05 and 0.05 for 
charged pions, kaons and protons respectively. The relative uncertainty due to stability 
evaluation of the fit is 0.06, 0.12 and 0.15 for charged pions, kaons and protons. 

Within the limited statistics of this sample the data are in good agreement with the 
generator prediction as well as with the prediction from MLLA (equation (4)). 

7.2 Average multiplicity 

After the event selection the analysis proceeds along the same lines as described for the 
qq analysis. Fully corrected ^-distributions are obtained and afterwards integrated to 
obtain the multiplicity. The results for the average multiplicity are shown in table 5 and 
compared to the predictions from PYTHIA 5.7 and HERWIG 5.8, with the systematic 
errors estimated as in subsection 6.7. The relative size of the systematic uncertainty 
due to extrapolation into the unseen region is 0.023 (0.024), 0.0023 (0.0022) and 0.0011 
(0.0011) for charged pions, kaons and protons from fully hadronic (semi-leptonic) W 
decays. The relative size of the systematic uncertainty due to particle identification is 
0.0052 (0.0093), 0.026 (0.076) and 0.088 (0.59) for charged pions, kaons and protons 
from fully hadronic (semi-leptonic) W decays. Finally, the relative size of the systematic 
uncertainty due to uncertainty in the estimate of the subtracted background is 0.0013 
(0.0003), 0.0027 (0.0004) and 0.0022 (0.0011) for charged pions, kaons and protons from 
fully hadronic (semi-leptonic) W decays. 
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event 




Particle 


(n> 


type 






(Data) 


PYTHIA 5.7 


HERWIG 5.8 


WW - 


■> qqqq 




31.65 ±0.48 ±0.76 


31.74 


31.55 








4.38 ±0.42 ±0.12 


4.05 


3.63 






PP 


1.82 ±0.29 ±0.16 


1.89 


1.66 


WW - 


-> qqli/ 




15.51 ±0.38 ±0.40 


15.90 


16.08 








2.23 ±0.32 ±0.17 


2.02 


1.82 






PP 


0.94 ±0.23 ±0.06 


0.95 


0.83 



Table 5: Average multiplicity for 7r ± , K* 1 and protons for WW events at 189 GeV. In the 
data the first uncertainty is statistical and the second is systematic. 

The ratios of the average multiplicities in (4q) events to twice the multiplicities in (2q) 
events for different hadron species, conservatively assuming uncorrelated errors, for the 
full momentum range and for momentum between 0.2 and 1.25 GeV/c are shown in ta- 
ble 6, and are compatible with unity. The systematic errors for the restricted momentum 
range were assumed to be proportional to the sum in quadrature of the systematic errors 
in the full momentum range, excluding the contributions from the extrapolation. 



Particle 


(n)(4g)/(2»(2g)) 


Allp 


0.2 GeV/c <p< 1.25 GeV/c 




1.02 ±0.03 ±0.04 


1.03 ±0.03 ±0.01 




0.98 ±0.17 ±0.08 


0.96 ±0.38 ±0.08 


pp 


0.97 ±0.28 ±0.11 


0.72 ±0.57 ±0.08 



Table 6: Ratio of average multiplicities in (4q) events to twice the values in (2q) events 
for 71^, and protons from WW events at 189 GeV, for different momentum ranges. 
The first uncertainty is statistical and the second is systematic. 



8 Summary and Discussion 

The mean charged particle multiplicities (n) and dispersions D in qq events at the 



different centre-of-mass energies were measured to be: 

(rc)i83GeV = 27.05 ± 0.27(stat) ±0.32(syst) (29) 

^issgcV = 8.08 ±0.19(stat) ±0.U(syst) (30) 

(n)i89GeV = 27.47 ± 0.18(stat) ±0.30{syst) (31) 

£>i89GeV = 8.77±0.13(staf)±0.11(sj/st) (32) 

(™>2oogcV = 27.58 ± 0.l9(stat) ±0A5(syst) (33) 

Afflocev = 8.64±0.13(stat) ±0.20{syst). (34) 

Figure 15 shows the value of the average charged particle multiplicity in e + e _ — > qq 



events at 183, 189 and 200 GeV compared with lower energy points from JADE [39], 
PLUTO [40], MARK II [41], TASSO [42], HRS [43], and AMY [44], with DELPHI results 
in qq7 events at the Z° [45], with the world average at the Z° [29], and with LEP results 
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at high energy [46-54]. The points from JADE, PLUTO and MARK II do not include 
the decay products of short lived K° and A. The value at M z o has been lowered by 0.20 
to account for the different proportion of bb and cc events at the Z° with respect to the 
e + e~ — > 7* — > qq [55]. Similarly, the values at 133, 161, 172, 183 and 189 GeV were 
lowered by 0.15, 0.12, 0.11, 0.08, 0.05 and 0.07 respectively. The QCD prediction for 
the charged particle multiplicity has been computed as a function of a s including the 
resummation of leading (LLA) and next-to-leading (NLLA) corrections [56]: 



(35) 



where s is the squared centre-of-mass energy and a is a parameter (not calculable from 
perturbation theory) whose value has been fitted from the data. The constants b = 0.49 
and c = 2.27 are predicted by theory [56] and a s (-\/s) is the strong coupling constant. 
The fitted curve to the data between 14 GeV and 200 GeV, excluding the results from 
JADE, PLUTO and MARK II, is plotted in Figure 15, corresponding to a = 0.045 and 
<^ s ( m z) = 0.112. The multiplicity values are consistent with the QCD prediction on the 
multiplicity evolution with centre-of-mass energy. 

The ratios of the average multiplicity to the dispersion measured at 183 GeV, 189 GeV 
and 200 GeV, (n) / D = 3.35 ±0.11, (n)/D = 3.14 ±0.07 and (n)/D = 3.19 ±0.10 (where 
the errors are the sum in quadrature of the statistical and of the systematic) respectively, 
are consistent with the weighted average from the measurements at lower centre-of-mass 
energies (3.13 ± 0.04), as can be seen in Figure 16. From Koba-Nielsen-Olesen scaling [57] 
this ratio is predicted to be energy- independent. The ratio measured is also consistent 
with the predictions of QCD including 1-loop Higher Order terms (H.O.) [58]. 

For WW events the measured multiplicities in the fully hadronic channel are: 

(n (4<?) )i89GcV = 39.12 ± 0.33(stat) ±0.36(syst) (36) 
(™ (4<?) )i83GeV = 38.11 ±0.57(stat) ±0.U(syst), (37) 

while for the semileptonic channel they are: 

(n (2<?) )i89GeV = 19.49 ± 0.31(stat) ±0.27(syst) (38) 
(n (2q) )i83GeV = 19.78 ± 0A9(stat) ± 0A3(syst) . (39) 

The PYTHIA Monte Carlo program with parameters tuned to the DELPHI data at 
LEP 1, predicts multiplicities of 38.2 and 19.1 for the fully hadronic and semileptonic 
events respectively. 

A possible depletion of the multiplicity in fully hadronic WW events with respect to 
twice the semileptonic events, as predicted by most colour reconnection models, is not 
observed in the full momentum range, in agreement with the results from other LEP 
collaborations [59]: 

/ n (4<?)\ 

* = 0.990 ±0M5(stat) ±0Ml(syst), (40) 



nor in the momentum range 0.1 < p < 1.0 GeV/c: 

( n (4q)\ 0.1<p<lGeV/c 

^P^yy = 0.981 ±0.024(stat) ±0M3(syst). (41) 
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No significant difference is observed between the dispersion in fully hadronic events 
and \/2 times the dispersion in semileptonic events: 

£)(4g) 

^ D(2q) = 0.94±0.03(stat)±0.03(syst). 



Assuming uncorrelated systematic errors for the two centre-of-mass energies, 183 and 
189 GeV, the inverse of the sums in quadrature of the statistical and the systematic errors 
were used as weights to give the averages: 

( n ( w )) = 19.44 ± O.lS(stat) ±0.12(syst) (42) 
D (W) = 6.20±0.11(stat)±0.06(sj/s*)> (43) 

where the multiplicities and their errors in (4q) were divided by 2 and the dispersions 
and their errors in (4q) were divided by \/2. 

The value of (rS w ^) is plotted in Figure 15 at an energy value corresponding to the W 
mass, with an increase of 0.35 applied to account for the different proportion of events 
with a b or a c quark. The measurement lies on the same curve as the neutral current 
data. The value of (n^) / = 3.14 ± 0.07(stat + syst), plotted in Figure 16, is also 
consistent with the e + e~ — > 7* — > qq average. 

The production of ir + , K + , K°, p and A from qq and WW events at 189 GeV was also 
studied. The results on the average multiplicity of identified particles and on the position 

of the maximum of the £ p = — log(^=) distribution were compared with predictions of 
PYTHIA and with calculations based on MLLA+LPHD approximations. Within their 
uncertainties the data are in good agreement with the prediction from the generator as 
well as with the predictions based on the analytical calculations in the MLLA framework. 

The ratio of the multiplicities of identified heavy hadrons in (4q) events to twice those 
in (2q) events is compatible with unity, both for the full momentum range and momenta 
between 0.2 and 1.25 GeV/c: 

Particle All p 0.2 GeV/c < p < 1.25 GeV/c 

K± 0.98 ±0.17 ±0.08 0.96 ± 0.38 ± 0.08 
pp 0.97 ±0.28 ±0.11 0.72 ±0.57 ±0.08 . 

The evolution for the £* for identified hadrons is in good agreement with the prediction 
from perturbative QCD (equation (4)). This underlines the applicability of MLLA/LPHD 
for the description of hadron production in e + e~ annihilation over the full LEP energy 
range. 
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Figure 1: Charged particle multiplicity distributions for (a) the (4q) events and (b) the 
(2q) events at 183 GeV, for (c) the (4q) events and (d) the (2q) events at 189 GeV. 
The error bars in the data represent the statistical errors. The shaded areas represent 
the background contribution; the histograms are the sum of the expected signal and 
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Figure 2: (a) Corrected momentum distributions of charged particles for (4q) events 
(closed circles) and (2q) events (open circles), compared to simulation without colour 
reconnection, at 189 GeV. The error bars in the data represent the statistical errors. The 
difference between (4g) and twice (2g) is shown in (b). 
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Figure 3: (a) Corrected momentum distributions of charged particles for (4g) events 
(closed circles) and (2g) events (open circles), compared to simulation without colour 
reconnection, at 183 GeV. The difference between (4g) and twice (2q) is shown in (b). 
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Figure 4: (a) Corrected £e distributions of charged particles for (4q) events (closed cir- 
cles) and (2q) events (open circles), at 189 GeV, compared to simulation without colour 
reconnection. The error bars in the data represent the statistical errors. The difference 
between (4g) and twice (2q) is shown in (b). 
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Figure 5: (a) Corrected £e distributions of charged particles for (4q) events (closed cir- 
cles) and (2q) events (open circles), at 183 GeV, compared to simulation without colour 
reconnection. The error bars in the data represent the statistical errors. The difference 
between (4g) and twice (2q) is shown in (b). 



32 



> 

CD 

o 

o 

H 
Oh 



4=J 
o 



> 

o 

in 
o 

H 
Oh 



4=1 
o 



3 -q- 

2.5 -E 

2 -E 

1 .5 -= 

1 -E 

0.5 -2 

-2- 

-0.5 -E" 

-1 -2 

-1 .5 -| 

-2 -2 

-2.5 -2 

-3 -- 




DELPHI 

1 89 GeV 




p T (GeV/c) 



DELPHI 

1 89 GeV 
b) 




• WW - 2W 
□ PYTHIA 



4 



p T (GeV/c) 



Figure 6: (a) Corrected px distributions of charged particles for (4g) events (closed circles) 
and (2q) events (open circles), compared to simulation without colour reconnection, at 189 
GeV. The internal error bars in the data represent the statistical error and the external 
error bars represent the sum in quadrature of the statistical and systematic errors. The 
difference between (4g) and twice (2q) is shown in (b). 
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Figure 7: (a) Corrected px distributions of charged particles for (Aq) events (closed circles) 
and (2q) events (open circles), compared to simulation without colour reconnection, at 183 
GeV. The internal error bars in the data represent the statistical error and the external 
error bars represent the sum in quadrature of the statistical and systematic errors. The 
difference between (4g) and twice (2q) is shown in (b). 
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Figure 8: Fractions of identified particles in radiative Z° events from 189 GeV data. 
The data have been taken under the same conditions as the signal data. Data (points) 
are in good agreement with the prediction from JETSET (lines) including full detector 
simulation. The top line and points indicate the sum of the fractions after unfolding. 
The error bars represent only the statistical error. 
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DELPHI preliminary 





Figure 9: £ p distributions (efficiency corrected and background subtracted) for charged 
particles, pions, kaons and protons in qq events at 189 GeV. Data (points) are compared 
to the prediction from JETSET (solid line). Only the statistical uncertainties are shown. 
The dashed dotted line shows a fit to equation (3). 



30 




Figure 10: £ p distributions (efficiency corrected and background subtracted) for charged 
particles, pions, kaons and protons in fully hadronic WW events at 189 GeV. Data (points) 
which have been boosted back to the W rest-frame, compared to the prediction from 
JETSET (solid line). Only the statistical uncertainties are shown. The dashed dotted 
line shows a fit to equation (3). 
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Figure 11: £ p distributions (efficiency corrected and background subtracted) for charged 
particles, pions, kaons and protons in semileptonic WW events at 189 GeV. Data (points) 
which have been boosted back to the W rest-frame, compared to the prediction from 
JETSET (solid line). Only the statistical uncertainties are shown. The dashed dotted 
line shows a fit to equation (3). 
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Figure 12: £ p distributions (efficiency corrected and background subtracted) for neutral 
Kaons and neutral Lambdas in e + e~ — > qq at 183 GeV and 189 GeV respectively, for 
data (points), and simulation using JETSET (dashed-dotted histogram). The full curves 
show the fit of the data to a Gaussian and the dashed ones to a distorted Gaussian. The 
error bars represent only the statistical error. 
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Figure 13: Average multiplicity of K + (a), K° (b), p (c) and A (d) as function of the 
centre-of-mass energy. Black squares are DELPHI high energy data, open squares are 
measurements from a previous DELPHI publication [34], and open bullets are values 
taken from the PDG [29] . Simulations using PYTHIA tuned to DELPHI data (solid line) 
and HERWIG 5.8 (dashed line) are superimposed. The error bars represent the sum in 
quadrature of the statistical and the systematic uncertainties. 
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Figure 14: Evolution of £* in e + e _ — > qq with increasing centre-of-mass energy. This 
measurement is compared with previous measurements from DELPHI from references [34] 
and lower energy experiments [37]. The lines are fits to equation (4) for the hadrons listed 
in the figure. The error bars represent the sum in quadrature of the statistical and the 
systematic uncertainties. For better legibility the data points have been shifted from the 
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Figure 15: Measured average charged particle multiplicity in e + e~ — > qq events as a 
function of centre-of-mass energy ^fs. DELPHI high energy results are compared with 
other experimental results and with a fit to a prediction from QCD in Next to Leading 
Order. The error bars represent the sum in quadrature of the statistical and the sys- 
tematic uncertainties. Some points are slightly shifted on the abscissa for clarity. The 
average charged multiplicity in W decays is also shown at an energy corresponding to the 
W mass. The measurements have been corrected for the different proportions of bb and 
cc events at the various energies. 
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Figure 16: Ratio of the average charged particle multiplicity to the dispersion in 
e + e~ — > qq events at 183, 189 and 200 GeV, compared with lower energy measurements. 
The error bars represent the sum in quadrature of the statistical and the systematic un- 
certainties. Some points are slightly shifted on the abscissa for clarity. The ratio in W 
decays is also shown at an energy corresponding to the W mass. The straight solid and 
dotted lines represent the weighted average of the data points and its error. The dashed 
line represents the prediction from QCD (see text). 



